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Abstract 

A QCD based calculation of the photon spectrum in semi-inclusive r]Q 
decays is performed. The method applied is an effective theory based on 
a IjmQ expansion of QCD. 





1 Introduction 


The expansion of QCD in inverse powers of the mass has become an indispensable 
tool for the analysis of heavy hadron decays. This method, formulated as an 
effective held theory known as Heavy Quark Effective Theory (HQET) [|I|, has 
put the physics of heavy quarks on a QCD based and thus model independent 
framework. 

While applications of these methods to systems with a single heavy quark has 
been studied extensively (see 0 for reviews), there is a lot of motivation to extend 
this ideas to systems with two or more heavy quarks, because heavy quarkonia 
(the and the T’s) have been known since some time and the observation of 
havored “doubly heavy” systems (such as the family of states) is very likely 
in the near future. 

The starting point of the construction of an effective theory for these systems 
is again the expansion of the relevant Greens functions in l/mg, where rriQ is 
the mass of the heavy quark. However, it turns out [Q that one can not use the 
static limit for two heavy quarks, in particular if their velocities differ only by 
an amount of order l/mg, i.e. vv' — 1 ~ ^qcd/'^q as it is the case in a particle 
containing two heavy quarks. The static limit breaks down and one is forced to 
include at least the kinetic energy into the leading order dynamics. 

However, unlike the case of one heavy quark an effective theory for heavy 
quarkonia in terms of the quark and gluon degrees of freedom is more complicated. 
This is mainly due to the fact that in a heavy quarkonium not only the scale of the 
light degrees of freedom (as e.g. A in HQET) appears as a scale small compared 
to the heavy quark mass, but also scales which do depend on the heavy quark 
mass. This point becomes clear if one studies a system like the positronium first. 

In the QED case one may safely assume that the coupling a does not run 
when one studies a QED bound state consisting of two heavy fermions, which we 
shall assume to be of the same mass mq. The coupling a is a small parameter 
and hence the inverse Bohr radius mqa of such a system is small compared to 
the heavy mass mq. Furthermore, in the case of QED it makes sense to talk 
about the binding energy Eb of this system, which is again smaller by a factor 
a, Eb ~ mqo^. 

Thus in a coulombic system there are two scales, the inverse Bohr radius and 
the binding energy which are both small compared to the heavy mass: 

mqo^ < mqa < mq. (1) 

In addition, both scales depend on the heavy mass, in this simple example in 
a trivial way. The expansion parameter in this case is the relative velocity Vrei 
between the two heavy constituents, and a consistent scheme for the velocity 
counting has been set up. The velocity in a coulombic system turns out to be 
Vrei = and hence is indeed a small parameter to expand in. In particular, the 
binding energy is of the order 
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Switching now to full QCD the coupling starts to run, introducing another 
scale Aqcd- The problem in applying the coulombic hierarchy of scales to realistic 
quarkonia is related to the question at which point in the hierarchy ([^) Aqcd 
enters. This is obviously a question how large thq is, and for “superheavy” 
quarkonia one can have the scenario Aqcd "C rriQa^ such that the binding is 
coulombic and - from the point of view of QCD - is of perturbative nature, 
since agiEs) <C 1. This case is governed by Non-Relativistic QCD (NRQCD) as 
formulated by Bodwin, Braaten and Lepage [Q, where the velocity counting of 
the QED analogue is used. In fact, this type of counting is also correct as long as 
the binding is due to instantaneous gluons such that the emission of a dynamical, 
transverse gluon is suppressed by some power of the velocity. 

Realistic quarkonia, such as charmonium of bottomonium, do not seem to 
fulhll the requirement of being coulombic. This becomes obvious from the spectra 
as well as from the leptonic widths, which are related to the wave functions at 
the origin in the simplest picture. Here it seems that Aqcd is at least larger 
than the binding energy or maybe even larger than the inverse Bohr radius. In 
such a case it is not obvious what the parameter is in which an expansion should 
be formulated. The I/tuq expansion yields in the numerators combinations of 
covariant derivatives, some of which may be identihed with the relative velocity. 
In a “minimalistic” approach (called HQQET in [^]), all derivatives are kept at the 
price of a proliferation of unknown matrix elements. In other words, in HQQET 
simply the powers of the operators are counted and one can translate NRQCD 
into HQQET and vice versa. Ending that some matrix elements in HQQET at 
leading order are down by factors of Vrei in the counting scheme of NRQCD. 

In the present paper we shall study a simple process which could shed some 
light on the question of how to correctly count the “powers” of a matrix element 
in the l/mg expansion for realistic heavy quarkonia. We shall discuss the decay 
7 + light hadrons, in which the hard subprocess is the annihilation of 
two heavy quarks into a photon and a gluon, which implies that the two heavy 
quarks have to be in a color octet state at the matching scale tuq. The matching 
to the eEective theory thus yields (at the matching scale fi = uiq) four-quark 
operators in which the quark-antiquark pair is in a color octet state, and hence 
in naive factorization these would vanish. Also in NRQCD these matrix elements 
would require at least an additional dynamical gluon, the emission of which is 
suppressed by factors of Vreu hence this process is an interesting laboratory to test 
the size of these octet contributions which seem to play some role in quarkonia 
production. 

In the next section we give a brief review of the efiective theory approach used 
here. In section 3 we set up the necessary machinery for these decays in which 
the two heavy quarks annihilate and review the short distance contributions. In 
section 4 we discuss the hadronic matrix elements and the restrictions imposed by 
heavy quarkonia spin symmetry. In section 5 we apply this method to calculate 
the total rate and the photon spectrum in 7 -|- light hadrons. Finally we 
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discuss our results and conclude. 


2 Structure of the Effective Theory for Quarko- 
nia 


As we shall exploit the fact that the mass of the heavy quark is large compared 
to both Aqcd and all other scales such as the analoga of the inverse Bohr radius 
and the binding energy, we are starting from a l/mg expansion of the QCD 
Lagrangian and the corresponding expansion of the fields up to order l/rriQ. It 
is well known that there are different ways to perform such an expansion (as e.g. 
by integrating out the small component fields from the QCD functional [§] or by 
a Foldy-Wouthuysen Transformation |^); although the I/tuq expansions of their 
Lagrangians and of their helds look different, the result for the Greens functions 
will be the same. Hence one may pick the most convenient representation for the 
application under discussion and we pick 
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where v is the velocity of the quarkonia, and the transverse components of the 
derivative is given by — VgVy)D^. The Ki and Mj are operators of 
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Note that at order l/mg one has also a purely gluonic contribution due to closed 
heavy quark loops. It contains the gluonic field strength tensor which is 
defined by igsG^^ = [Dfi, Dt,]. 

The first terms of these two expansions (@) define the static limit, which 
has been successfully applied to systems with a single heavy quark. In order to 
describe a system with more than one heavy (anti)quark one has to write down 
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the same expansion (2) for each heavy quark. However, there is no static limit 
for a system with two heavy quarks if the two heavy quarks move with almost 
the same velocity as it is the case for a quarkonium; one runs into problems with 
diverging phases and “complex anomalous dimensions”, which are considered in 
detail in [Q. 

In order to cure this problem one has to choose the unperturbed system such 
that these phases are already generated by the leading order dynamics, i.e. instead 
of the static limit one has to use the non-relativistic Lagrangian. In this case the 
“power counting” has to be modihed in such a way that the leading term is the 
static plus the kinetic energy term, leading to the fact that the “time derivative” 
ivD has to be counted as two powers of the “spatial derivative” iD^. In physical 
terms, for a bound state this corresponds to the balance between potential and 
kinetic energy, since the covariant time derivative ivD contains also the potential. 

For a system with a heavy quark and a heavy antiquark one then starts from 

Co = hi+\ivD)hi^^ - hi-\ivD)hi-^ + (4) 

zniQ ^niQ 

where we have assumed for simplicity that the two quarks have the same mass; 
the case of unequal mass is obvious. 

Naive power counting would suggest to include also the Pauli term a ■ B into 
the leading-order Lagrangian. However, in a coulombic system this term is in 
fact down by powers of Vrei, since the chromomagnetic held is generated by the 
relative motion of the two heavy quarks; this has been made manifest in [§]. 

Most of the success of HQET is due to heavy quark havor and spin symmetry. 
However, once one uses (^) the symmetries are somewhat different for HQQET or 
NRQCD. First of all, (|[) depends on the mass through the kinetic energy term; 
consequently the states will depend on rriQ in a non-perturbative way and heavy 
havor symmetry is lost. In this way the mass dependent small scales (the inverse 
Bohr radius and the binding energy) are generated by taking matrix elements 
of operators involving derivatives with these states. On the other hand, since 
the chromomagnetic moment operator does not appear to leading order, (^) does 
not depend on the spins of the two heavy quarks, and hence there is a spin 
symmetry which is larger than in HQET because we have two heavy quark spins; 
the resulting symmetry is an SU{2) ® S'f/(2) corresponding to separate rotations 
of the two spins. 

For the case of heavy quarkonia all states fall into spin symmetry quartets 
which should be degenerate in the non-relativistic limit. In spectroscopic notation 
these quartets consist of the states 

[rCif, n^ie n^ie+i]- ( 5 ) 

For the ground states the spin symmetry quartet consists of the tjq (the 0“ state) 
and the three polarization directions of the Tq (the 1“ state). 
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The heavy quarkonia spin symmetry restricts the non-perturbative input to 
a calculation of processes involving heavy quarkonia. Of particular interest are 
decays in which the heavy quarks inside the heavy quarkonium annihilate. The 
annihilation is a short distance process that can be calculated perturbatively in 
terms of quarks and gluons, while the long distance contribution is encoded in 
certain matrix elements of quark operators. Logarithmic dependences on the 
heavy quark mass may be calculated by employing the usual renormalization 
group machinery. 


3 Annihilation Decays of Heavy Quarkonia 

The starting point to calculate processes like rjQ —>■ light hadrons, ?7q —>■ 7 + light 
hadrons, or the corresponding decays of the Tq states is the transition operator 
T for two heavy quarks which annihilate into light degrees of freedom. This will 
in general be bilinear in the heavy quark fields, such that 

T{x, 0 = (-^)g(x + OK{x, OQ{x - 0 , (6) 

where K {X, ^) involves only light degrees of freedom and X and ^ correspond to 
the cms and relative coordinate respectively. If we identify the field Q with the 
quark and Q with the antiquark, so we shall make the large scale tuq explicit by 
redefining the fields as 

Q{x) = e'x.Y>{-imQ{vx))Q[l^\x), Q{x) = eyrp{-imQ{vx))Q[^\x) . (7) 

This corresponds to the usual splitting of the heavy quark momentum into a large 
part ruQV and a residual piece k. Inserting this into one finds 

T= (-^)expH2mQ^;X]g(-)(X + 0i^(X,0gi+HX-0. (8) 

The inclusive rate for the decay of a quarkonium T —> light degrees of freedom is 
then given by 


T 


(Tl 


d"XdWT(X,0Tt(0,O+h.c. 


IT). 


(9) 


The next step is to perform an operator product expansion for the non-local prod¬ 
uct of the quark field operators. This expansion will yield four-quark operators of 
increasing dimension starting with dim-6 operators. The increasing dimension of 
these operators will be compensated by inverse powers of the heavy quark mass, 
so generically the rate takes the form 

r = moEf-^') c(o7.a)('I'|o 7|'I'>U. (lo) 

n,i \^Q/ 
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where n = 6, 7,... is the dimension of the operator and i labels different operators 
with the same dimension. The coefficients are related to the short 

distance annihilation process and hence may be calculated in perturbation theory 
in terms of quarks and gluons. Once QCD radiative corrections are included, the 
IJ,) acquire a dependence on the renormalization scale /i which is governed 
by the renormalization group of the effective theory. The rate T is independent 
of and hence the /r dependence of has to be compensated by a 

corresponding dependence of the matrix elements. 

The non-perturbative contributions are encoded in the matrix elements of 
the local four-quark operators, and the mass dependence of these operators is 
also expanded in powers of l/mg and thus the remaining mq dependence of the 
matrix elements is only due to the states. In terms of the uiq independent static 
fields 


= h[^\x) + 0{l/mQ) , Q[ = h[ \x) + 0{l/mQ) ( 11 ) 

one has in total four dim-6 operators 

and 

with the color matrix C, where one has the two possibilities C®C = 1 ® 1 (color 
singlet) oi C ® C = ® T°‘ (color octet). These operators do not mix under 

renormalization, all anomalous dimensions vanish. 

There are no dim-7 operators, since these are either proportional to {ivD) 
and can be rewritten in terms of dim-8 operators by the equations of motion (see 
(^) or their forward matrix elements are forbidden by symmetries. At dim-8 one 
finds 30 local operators 

bP = [iDPhpPChP}] [iDphi-pchP)] , 
bP = itophppchp)] ItDlppPChpl)] , 

Bf> = liD^ChPllChP)] liDphPPChP)], 

cP = [iDl(hPl,ChP)] p-P(WpChP] + h.c,, 
cP = [iD'phPpChP] |A<->(i77)CA<+>| + iLC., 

cP = [iD'phpp^ChP)] PPP(WpChP] + hx., 

Cf > = \iDpUpiiChp^)\ VphpiDpChP] + hx., 

dP = WPiuPqchP] \hp^,(wpchfp] , 
dP = [hP{ip^)ChP] [hP{iip^)ChP\, 
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dP = [hi*hl(zD^)Chi-'l IIii->7piD;!r)ChP] , 

Cf > = lhi->7Hiip)ChP] , 

-Bp’ = \li[-l 7 ^{iD^YCh‘-* 1 \ + h.c.. 

Bp’ = [hPllChP] [Sp>(»irp(»B’i-)Cfcl+>] + h.c,, 

Bp’ = m*hlChl->] [S<-)(»B’i:)(»iP)Cfcp'] + h-c,, 

Bp’ = IfcPPiCfc!,"’] [Sppp (!-D-^)^Cfc!,+>] + h.c.. 
where we have defined 

f{x){D)g{x) = f{x){Dg{x)) - {Df{x))g{x) (12) 

corresponding to (the gange invariant generalization of) the relative momentnm. 

Note that we have split these operators into fonr blocks. The Bl ^ operators 
contain two derivatives acting on both the heavy qnark and the heavy antiquark 
operators; thus they correspond to a motion of the cms-system of the two heavy 

(C) 

quarks relative to the cms-system of the quarkonium. The operators of the Cl 
type are chosen such that one derivative acts again on both fields, while the other 
derivative corresponds to the relative motion of the heavy quark-antiquark pair. 
Finally, the and correspond to two powers of the relative momentum 
of the two heavy quarks. 


4 Renormalization Group Evolution 


In (^^ a spurious dependence on the renormalization point /i appears, where 
/i separates long from short distances. The coefficients in front of the matrix 
elements may be calculated in perturbation theory and are determined at the 
scale fi = rriQ from a matching calculation. Once the matching is done one has to 
evolve down to a small scale which is typical for the matrix elements appearing 
in ( 0 ). As it was discussed in the introduction there are several small scales in a 
quarkonium and in this context several questions arise. Firstly, if the dimension 
is not the correct way of counting the relevance of a matrix element, how is the 
counting of powers done correctly? Secondly, how far should one evolve down 
using the renormalization group in order to get the “best” result? 

For superheavy quarkonia, i.e. the ones one may obviously describe using 
NRQCD, these problems are easier to discuss. Since the binding is due to in¬ 
stantaneous coulomb gluons, one may write the Lagrangian in Coulomb gauge 
as 


C = hi+\tvD)hi+^ - hi-\ivD)hi-^ + 


2m, 
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2m, 


(13) 
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since the gluon fields associated with the spatial derivatives iD^ are suppressed 
due to the smallness of the coupling agirnQvl^i). Furthermore, coulomb gluons do 
not appear as dynamical degrees of freedom, rendering the superheavy quarko- 
nium as a real non-relativistic two particle state, with a potential determined from 
the gluon Greens function ( 0 |y 4 o(a;) 74 o( 2 /)| 0 ). In such a case one may perform a 
Fock state decomposition of the quarkonium of the form 

1 ^) = Co\QQ) + Ci\QQg) + C 2 \QQgg) + ■ ■ ■ (14) 


in which subsequent coefficients are suppressed by powers of Vreh Ci since 

each dynamical gluon is suppressed by a power of Vrei- In addition, one has the 
same selection rules as in Non-Relativistic QED (NRQED), yielding a consistent 
scheme for an expansion in Vrei- 

iC) 

In this way to count powers the matrix elements of the operators B\ ^ and 

(O') 

Cl are down by additional powers of Vrei compared to their dimension, while the 
operators of the and type are of order {mqVrei)^ and thus their matrix 
elements constitute the first subleading corrections in the power counting scheme 
along the lines of NRQED. The same argument applies in for the operators in 
which the two quarks are in a color octet state. Since the quarkonium has to 
form a color singlet, these contributions have to be of the form of a higher Fock 
component and are thus suppressed by additional powers of Vrei- 

However, realistic quarkonia seem not to be close to this limit, and it is well 
possible, that the Vrei counting scheme fails for these systems. In such a case it 
seems to be useful to count the operators according to their dimension given by 
the powers of iD-^. This means that all the matrix elements of the operators 
Bl \ ...,El ^ will be taken as the subleading terms. In particular this means 
that the states of the effective theory are not simply non-relativistic two particle 
states in which the two heavy quarks are bound by a potential generated by 
instantaneous, coulomb type gluons. Rather they have to contain dynamical 
gluons as well such that all the Fock components in (pH]) are of comparable size. 
The physical picture is illustrated in fig.|^. 

The operators bI^\ ..., E^^'^ are all local operators corresponding to the I/ttiq 
expansion of the hard kernel. In addition to these local contributions one has also 
non-local terms originating from single insertions of the Lagrangian of order l/rriQ 
and from double insertions of the Lagrangian of order IjmQ. These contributions 
correspond to the corrections to the states of the effective theory 

Under renormalization the local dim -8 operators do not mix: only the double 
insertion of the kinetic energy operator of order l/mg {Ki = + k[ mixes 

into some of the above operators. Denoting this contribution as 


rjn{C) _ 


fd‘‘xd*yTlAf* Ki(x)Ki(y)] 


( 15 ) 


one obtains in one-loop renormalization group improved perturbation theory two 
sets of equations for the coefficients of the operators with the spin structure 75 ( 8)75 
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Figure 1: Illustration of the various Fock components. In the standard 

NRQCD picture the left hgure represents the leading contribution. 

ii 

= C(DP.mo) + f 

C(E?,^) = C(£P,mo) - 

C{Bf\^L) = C(sf',mo) - 33 ^"*2 cur. mg) In,,, 

= C(Df'.mg) + „/^„ C(r<‘>,mg)lii,,, (16) 

+ —-^- C{T[^\mQ)\nr], 

3 33 - 2nf ^ ^ ’ 

C(£r,n) = C(£f>,mg) - Y 33_^2„. *^<h”,'"Q)ln'). 

where n/ is the number of active flavors and r] = /cxs{rnQ)). Furthermore, 

the coefficients C{Tl;^\mQ) are the same as the ones for the dim -6 operators 
C{Ai yTUg) since the kinetic energy operator is not renormalized. 

The second set of equations is for the operators with spin structure 7 ^ ( 8 ) 7 ^ 
and due to heavy quarkonia spin symmetry one obtains the same equations; all 
other renormalization group equations are trivial. 

The calculation of the coefficients is based on perturbation theory, and hence 
as has to be sufficiently small. Although the hnal result does not depend on /i, 
there is still the practical question of how far one should run using the renormal¬ 
ization group. Clearly /i should be chosen to be one of the three, the inverse Bohr 
radius, the binding energy or Aqcd- In the worst case Aqcd is of the order of the 
inverse Bohr radius, and then the renormalization group evolution has to stop 
there. In the superheavy case, Aqcd is small compared to the binding energy and 
one may run down to the scale set by the binding energy, below this scale one 
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would need to switch to an effective theory involving hadronic instead of QCD 
degrees of freedom. 


5 Hadronic Matrix Elements and Spin Symme¬ 
try 

A calculation of an annihilation decay then involves to calculate the C{Ol^\n) 
at the scale fi = mq by matching the effective theory to full QCD. Once this 
is done, one may run down to some small scale /i. In the superheavy case fi is 
of the order of the binding energy of the heavy quarkonium, thereby resumming 
the well known logarithms of the form ln{mQ/fi) that appear in the calculations 
of decay rates of heavy p-wave quarkonia. As an example we shall study the 
decay Pg —^ 7 + light hadrons in the next section and consider the first nontrivial 
corrections to this mode. 

The matrix elements of the operators Bf"\ as well as the non-local 

terms are non-perturbative qnantities, which are constrained by heavy qnarko- 
nia spin symmetry. In order to exploit this symmetry, one may use the usual 
representation matrices for the spin singlet and spin triplet qnarkonia 

H^iy) = y/MP+-f5 for A = 0 , H^iv) = ^MP+e for ^ = 1 (17) 

where M ^ 2mQ is the mass of the heavy quarkonium and P+ = (l-|-^)/2. Using 
this one finds for the matrix elements of the dim -6 operators 

= a^^\n,^)G with G = Pi(H 2 s+iT)PTiT'H^s+i) 

Thus for each n and ^ and for each color combination one finds a single parameter 
for both the spin singlet and spin triplet quarkonium. 

Correspondingly for the dim -8 operators we get 

(vl/|[^D;[(hWrCh(-))] [iD^{hi-^rGh[^^W) = - v,v,)G 

(vl/|[hWr(^5^)Ch(-)] +h.c.|T) = c^^\n,^){g,,-v,v,)G 

{^0+^T{Wi)Ghi-^] [h^-^T\wl)Gh^+^] = - v^v,)G 

(vl/|[hWr(*3^)(^3^^)C'h(-)] [h^^T'ChW] +h.c.|T) = e^^\n,i)ig^,-v^v,)G 

For fixed valnes of n and ^ one finds that eight parameters are needed to describe 
the matrix elements of the dim -8 operators. 

These matrix elements are non-perturbative, but from vacuum insertion one 
is lead to assume 

0) = ■^|i?„o(0)P 3 > £) for £ 7 ^ 0 , 

87r 
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Figure 2: The hard subprocesses for > 7 + light hadrons. 

= ^\KM? » dW{n,() for 1, (18) 

e<‘)(n.O) = -ReK„(0)K,(0)] » el‘>(n, f) for f #0, 

TT 

(n, 0) S> (n, £) for all u, £, 

1) 3> (n, i) for all u, £, (19) 

(u, 0) S> (u, i) for all n, i , 

where Rni{r) is the radial wave function of the quarkonium. The same reasoning 
yields the expectation that (n, i) and (n, i) are small compared to the 
coefficients that are non-vanishing in vacuum insertion. 

In fact, the expectations from vacuum insertion are strengthened by the su¬ 
perheavy case in which vacuum insertion is true to leading order in the Vrei 

iC) 

expansion. Furthermore, the color octet as well as the operators of the 5) and 
Cl type are also down by powers of Vrei- 


6 The Photon Spectrum in t/q ^ 7 + light hadrons 


The short distance part of the process under consideration is given by the am¬ 
plitude obtained from the Feynman diagrams depicted in hg.H. The amplitudes 
are evaluated in QCD and then matched to HQQET by performing an expan¬ 
sion in l/mq. At the matching scale the spectrum, including the hrst non-trivial 
corrections, has the form 


dV 

dx 

with 


= g 


5(1 - x)A + {5(1 - x)Bo + 5'(1 - x)B, + 5"(1 - x)B 2 } 


( 20 ) 


G = 


S'Kas{mQ)aemQQ 


m 


Q 
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( 21 ) 












where x := E^/mq is the scaled photon energy, aem = 1/137 is the electromag¬ 
netic coupling and Qq is the charge of the heavy quark. At the matching scale 
/i = mq the coefficients of the ^-functions in (pOD are given in terms of the forward 
matrix elements discussed in the last section and non-local terms originating from 
the correction in the Lagrangian. One hnds from the matching calculation 


^ = {'nQ\^i\'nQ)\mQ ( 22 ) 

+ (-i) J d^x {r]q\T{{K2{x) + Cgiue{x))Af^}\r]q)\mQ (23) 

iVQ\T{Mi{x)Mi{y)Af'^}\T]q)\,nQ 

B, = l{vq\B[^^\vq)\m^ - (24) 

132 = ^{VQ\B^^^\VQ)\mQ (25) 


The renormalization group may now be used to run down to some small scale 
/i ~ AqcD, thereby extracting logarithms of the form log{mq/fi) from the matrix 
elements. Using the one loop result as given in the last section, one has at the 
small scale 


A 

Bo 


Bi 

B2 




. 56 

4- 

33 — 2nf 

128 1 


ln?7 


{r]Q\E[^^\r]q)\g - 


7- 


Inr] 


\nr]{7]Q\D\^^\r]q)\f,+ 


3 33 -2nf 33-2nf 

+ {-i) [ d'^x {riq\T{{K2{x) + Cgiue{x))Af^}\r]q)\ 


33 — 2nf 

\nT]{T]Q\D^^'^\T]Q)\ 


(26) 

{VQ\B[^'’\Vq)\g 


+ 




d^x d^y {yq\T{{M,{x)M,{y)Af^}\yq) \ 


^{Vq\B?\Vq)\i^ -\{VQ\E?\VQ)\f^ 


:{VQ\Bl^\VQ)\^ 


(27) 

(28) 
(29) 


Note that the renormalization group flow has induced two operators that have not 
been present at the matching scale. These type of logarithms has been observed 
already in the calculation of p-wave quarkonia some time ago |^; these have been 
hxed order calculation and the logarithms appear here as infrared singularities. 
In the effective theory approach they are generated by the renormalization group 
flow which in addition even resums these terms, since 


- Inr^ = In (l -h a^(mQ)^ In — 

V 27r mq^ 
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= a, mQ — In-+ - 

ZTT rriQ 2 


/ X ^0 , /i 

(^s{mQ)—\n - 

271 ruQ^ 


+ ... , 


( 30 ) 


where (3q = 11 — 2nf. 

The expression obtained for the photon spectrum contains ^-function and its 
derivatives which is obviously unphysicah The origin of this singular behavior is 
of kinematical nature: At the partonic level, the initial state quarks move with 
the same velocity and hence act like a single particle, which then decays into two 
massless objects. Thus this is a two particle decay and so the energies of the hnal 
state particles are hxed. 

In order to compare with the observed hadron spectrum one has to apply 
some “smearing” in the sense of JOT; in particular, if one calculates moments of 
the spectrum ([^) one obtains a sensible answer which may be compared to the 
observed spectrumQ. In other words the l/mq expansion yields an expansion of 
the spectrum in terms of singular functions 


rir °° 1 

:r = S - x) , 

dx n\ 

where the coefficients of the expansion are the moments 

-•i . dV 


(31) 


Mn = / dx—(l — x)'^. 
Jo dx 


(32) 


The zeroth moment is simply the total rate 


r = ^[^ + r^^o], 


while the hrst and the second moments are entirely of order l/mq 

Ml = QBi , M2 = gB2. 


(33) 


(34) 


As expected, the hrst and the second moment are entirely proportional to color 
octet contributions, while the renormalization group how induces the color singlet 
operator In factorization as well as in the NRQCD case all these opera¬ 

tors are down by at least two powers of Vrei and hence this process would be 
very much suppressed compared to allowed processes such as J/T —>• 7 -|- light 
hadrons. Hence this process is a nice testing ground for factorization and the 
power counting scheme of NRQCD. 

^ In fact, the situation is completely the same as in the endpoint of inclusive B —> or 

B Xsj decays ||^, where the l/mg expansion also only yields the moments of the observed 
spectrum. 
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7 Conclusions 


Unlike systems with a single heavy qnark quarkonia-like systems are much more 
difficult to describe. The simplicity of HQET is due to the fact that only a single 
small scale appears, which is independent of the heavy mass, and thus a static 
limit may be performed. In a quarkonium several small scales appear which 
depend in a non-perturbative way on the heavy mass. The superheavy case re¬ 
sembles very much a QED like system which is bound by coulombic forces. In 
such a system one has the inverse Bohr radius mqVrei and the binding energy 
mqvl^i as small scales, and if Aqcd ^ ^^e binding is in the sense per¬ 

turbative that it is due to one gluon exchange which takes in coulomb gauge the 
form of an instantaneous potential. 

In this superheavy case one may transfer practically all the knowledge of 
the description of positronium in NRQED to the QCD case. However, realistic 
quarkonia do not seem to be too close to this limit, since the charm and the 
bottom mass are too small to fulhll Aqcd ^ ^Q^'^rei- Thus the power counting 
scheme of NRQED cannot be naively transferred to QCD, and the way how to 
organize the effective theory calculation becomes obscure. 

The safe way in this case is to rely only on the dimension of the operators 
involved, thereby taking into account a proliferation of unknown parameters, 
which all would need to be determined from experiment. It is then of some 
interest to dehne observables which are sensitive to specihc matrix elements and 
hence may shed some light on what happens in realistic quarkonia. 

In the present paper we have studied the decay ^ 7 + light hadrons, which 
would be strongly suppressed compared to J/T j + light hadrons, if both J/T 
and rjq were superheavy. However, this process is not suppressed due to higher 
powers of l/rnq appearing in the heavy mass expansion, rather it is suppressed 
due to the quarkonia states, with which the matrix elements are taken. 

It turns out that the moments of the measured photon spectrum in rjq — 7 + 
light hadrons is sensitive to matrix elements which would be strongly suppressed 
for the superheavy case. Unfortunately, there are not yet data on this process, 
such that a check, whether it is really suppressed as predicted by the counting of 
powers of Vrei has to wait for future experiments. 


Acknowledgments 

We would like to thank Martin Beneke and Gerhard Schuler for discussions and 
Andreas Schenk for checking details of the calculations. This work was supported 
by Deutsche Forschungsgemeinschaft. 


References 


14 



[1] M. Voloshin and M. Shifman, Sov. J. Nucl. Phys. 45 (1987) 292 and 47 
(1988) 511; 

N. Isgur and M. Wise, Phys. Lett. B232 (1989) 113 and B237 (1990) 527; 
E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511; 

B. Grinstein, Nucl. Phys. B339 (1990) 253; 

H. Georgi, Phys. Lett. B240 (1990) 447; 

A. Falk, H. Georgi, B. Grinstein and M. Wise, Nucl. Phys. B343 (1990) 1. 

[2] H. Georgi: contribution to the Proceedings of TASI-91, by R.K. Ellis et al. 
(eds.) (World Scientihc, Singapore, 1991); 

B. Grinstein: contribution to High Energy Phenomenology, R. Huerta and 

M. A. Peres (eds.) (World Scientihc, Singapore, 1991); 

N. Isgur and M. Wise: contribution to Heavy Flavors, A. Buras and M. 
Lindner (eds.) (World Scientihc, Singapore, 1992); 

M. Neubert, SLAG-PUB 6263 (1993) (to appear in Phys. Rep.); 

T. Mannel, contribution to QCD-20 years later, P. Zerwas and H. Kastrup 
(eds.) (World Scientihc, Singapore, 1993). 

[3] B. Grinstein, W. Kilian, T. Mannel and M. Wise, Nucl. Phys. B363 (1991) 
19; 

W. Kilian, P. Manakos and T. Mannel, Phys. Rev. D48 (1993) 1321; 

W. Kilian, T. Mannel and T. Ohl, Phys. Lett. B304 (1993) 311. 

[4] G.T. Bodwin, E. Braaten and G.P. Lepage, Phys. Rev. D51 (1995) 1125. 

[5] T. Mannel and G. Schuler, Z. Phys. C67 (1995) 159; 

T. Mannel and G. Schuler, Phys. Lett. B349 (1995) 181. 

[6] T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. B368 (1992) 204. 

[7] J. Korner and G. Thompson, Phys. Lett. B264 (1991) 185; 

S. Balk, F. Korner and D. Pirjol, Nucl. Phys. B428 (1994) 499. 

[8] M. Luke and A. V. Manohar, UTPT 96-14, UGSD/PTH 96-24, |hep- 
I ph/9610^ 

[9] R. Barbieri, R. Gatto and E. Remiddi, Phys. Lett. B61 (1976) 465. 

[10] E.G. Poggio, H.R. Quinn and S. Weinberg, Phys. Rev. D13 (1976) 1958. 

[11] M. Neubert, Phys. Rev. D49 (1994) 3392; 

M. Neubert, Phys. Rev. D49 (1994) 4623; 

LI. Bigi, M.A. Shifman, N.G. Uraltsev and A.I. Vainshtein, Phys. Rev. D52 
(1995) 196. 


15 








